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ABSTRACT 


This  report  contains  the  development  of  the  theoretical 
equation  for  the  aerodynamic  torque  absorbed  by  a helicopter 
rotor  which  is  experiencing  transient  disturbances  in  RPM 
while  the  helicopter  is  hovering.  Inasmuch  as  the  dynamical 
equation  for  torque  was  found  to  include  the  effects  of  the 
simultaneous  transients  in  blade  coning  angle  and  vertical 
vnl-city  of  the  aircraft,  it  was  necessary  to  derive  auxiliary 
equations  defining  these  additional  variables,  and  these 
equations  are  also  shown  here. 

The  system  of  equations  derived  in  this  report  is  required 
in  the  analytical  study  of  controls  for  automatic  regulation 
of  helicopter  rotor  RPM. 
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I INTRODUCTION 


The  analytical  study  of  helicopter  rotor  RPM  variation 
and  the  automatic  control  thereof  is  based  to  a large  extent 
on  the  relationship  from  dynamics, 

- I*#  a (1) 

This  simply  states  that  the  unbalanced  torque  acting  on  the 
rotor  gives  rise  to  an  angular  acceleration  (about  the  rotor 
shaft)  whose  magnitude  also  depends  on  the  effective  moment  of 
inertia  of  -he  rotating  system.  The  value  of  the  aerodynamic 
torque  (sometimes  referred  to  as  the  ’required’  torque  or 
the  torque  ’absorbed' by  the  rotor)  has  been  determined  theoreti- 
cally for  steady-state  rotor  operation  by  many  investigators 
(e.g..  Reference  (1),  equations  2.29  and  4°H5)  » For  the  case 
of  variable  RPM,  however,  as  indicated  by  equation  (1),  the 
expression  for  has  apparently  not  been  previously  derived. 

The  derivation  of  this  expression  was  the  principal  objective 
of  the  study  reported  herein. 

The  conditions  stipulated  for  the  present  treatment  for 
the  determination  of  as  well  as  other  requirements  and 
assumptions  involved  in  this  work,  are  described  in  the  follow- 
ing: 

1.  A single,  main-rotor  type  helicopter  was  selected 
for  analysis. 

2.  It  is  presumed  that  the  helicopter  was  hovering 
at  a point  in  space,  and  that  it  subsequently 
experienced  a vertical  gust,  a change  in  collective 
pitch  (say  due  to  pilot  action),  ami/  or  a change 
in  rotor  RPM  due  to  an  engine  disturbance. 

3°  It  is  assumed*’*'  that  as  a consequence  of  the 
above  disturbances,  the  helicopter  moves  only 
in  a vertical  direction,  with  the  tip  path  plane 
remaining  horizontal. 

While  the  induced  velocity  at  the  rotor  was  held 
to  be  invariant  with  time  throughout  the  transient 
experienced  by  the  helicopter,  it  was  presumed 


**  This  will  be  shown  in  a later  work  to  be  a fact  rather  than 
an  assumptionp  at  least  when  the  tail  rotor  is  disregarded. 
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to  vary  with  blade  radius  in  the  manner  prescribed 
by  the  combined  blade  element-momentum  theory* 

This  is  in  accordance  with  the  conclusions  pre- 
sented in  Reference  (2).  Hence, 


Tip  loss  effects  were  evaluated  by  means  of  the 
well-known  factor  B in  computing  blade  lift  forces* 

5*  Aerodynamic  forces  were  assumed  to  develop  on 

the  rotor  blade  elements  without  time  lag*  These 
forces  were  also  presumed  to  be  independent  of 
Mach  pumber  and  Reynolds  number  variations  during 
the  transients* 

6.  To  avoid  encountering  non-linear  differential 
equations,  all  variations  of  parameters  were 
taken  to  be  of  very  small  magnitude* 

7*  Elastic  properties  of  the  rotor  were  ignored. 

8.  The  drag  coefficient  at  the  blade  element  was 
taken  to  be 

as  suggested  by  Sisslngh,  Bailey,  Qessow  (see, 
for  example.  Reference  (3),  page  8l),  and  others* 

9*  The  rotor  blades  were  presumed  rectangular  and 
untwisted,  and  the  flapping  hinge  offset  e was 
considered  negligibly  small.  hinge  and  drag 
hinge  effects  were  not  considered. 
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II  DERIVATION  OP  THE  TORQUE  EQUATION 


Following  the  usual  approach  for  the  determination  of 
the  aerodynamic  forces  on  the  rotor,  based  on  blade  element 
concepts,  the  diagram  of  the  velocities  and  elemental  forces 
is  considered*: 


9 o(t 


Up\± £ 


4-^fc. 


Figure  1 


By  examination  of  Figure  1,  it  is  evident  that  the  elemental 
aerodynamic  torque  about  the  rotor  axis  for  the  blade  element 
at  radius  r is  given  by 


c/Q/)  •rdL  sin  (p  +HO  ras  (() 


or,  since  Up  is  everywhere  much  smaller  than  Uip, 

o/Qt-r-fedL+rdD 


For  the  hovering  regime, 

/ / D _ **  / D 


Ur  = r£l  - r(±U  +D.&) 
Up  - ~ 4 


See  Section  VI  for  explanation  of  the  symbols  employed  here: 
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Note  that  equation  (7)  provides  for  the  possibility  of  a vertical 
gus t of  velocity  w being  imposed  on  the  rotor,  the  (+)  sign 
implies  that  .the  gust  is  directed  downwardly*  The  (-)  sign 
in  front  of  (which  is  the  vertical  transient  velocity  of  the 
helicopter)  accounts  for  the  convention  that  downward  movements 
of  the  aircraft  are  positive.  The ^disturbance  variables  In 
equation  (6)  and  (7)»  , p and  4.  , do  not  vary  with  blade 

radius,  and  are  small  quantities « ^The  gust  velocity  w which 
produces  the  transient  effects  is  also  presumed  to  be  small. 

Also  from  Figure  (1), 


o',  - 0-0 


(8) 

(9) 


It  is  presumed  that  the  lift  coefficient  is  given  by  the  usual 
expression 


Ct  = a(at) 


Consequently, 


dQj-  @f-UzCj_dr  +■  r ^UzC0odr 

- j I 


(10) 


(ID 


-Pt 


UpUj-aCO-tP)  +Up(So  +£,&r  *Sz0cf^rdro.z) 


Making  the  appropriate  substitutions,  neglecting  non-linear 
and  other  second  order  terms,  and  integrating  (to  the  tip 
radius  R in  the  case  of  drag  terms,  and  to  BR  in  the  case  of 
lift  terms) , 


0,(4,  * 6,/S  + C,  ur -6,  j-c)  +0/1, 

- 


€,  + 2/1, ft  +2D,ur-?Qjz 
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+e2  Cl1  \so  +S,(90  +92)  +S2  (Ol  +29e92) 


+3Z(£Q.0Q2)  §o+St00  + &z9q 


-AzEl0  S,  + zSz  (90  +92) 


-S,+2Sz0o  Cl0(3z^  + Czar-Cz^2)^Az2U 


+8Z  Ez  +2/izfi  +2Dzur  -ZDzfc  +<i&AQt>A,  f Od) 


where  the  following  symbols  have  been  used  to  represent  the 
separate  integrals  involved  in  obtaining  equation  (13): 

,*e 

A j r'lrdr 

6>-Ajpf*r'dr 

l r*dr 
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Uae  of  the  subscript  2 instead  of  1 with  the  letter  symbols  in 
equation  (14)-(19)»  as  for  example,  A£,  B2»  etc*,  means  that 
the  upper  integration  limit  is  R instead  of  BR. 

As  indicated  in  equation  (13),  the  aerodynamic  torque 
during  the  transient  depends  not  only  on  £2^  » the  time  variant 
disturbance  in  rotor  speed,  but  also  on  ^3  and  4C  , the 
blade  flapping  and  helicopter  vertical  motion,  respectively. 
Thus,  for  evaluation  of  Q.A,as  required  for  the  solution  of  equa 
tion  (1),  it  is  necessary  to  write  simultaneous  equations  de-r 
fining  the  flapping  and  helicopter  vertical  motion.  This  is 
the  subject  of  the  next  two  sections  of  this  report. 
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Ill  THE  BLADE  CONINQ  EQUATION 


Consider  the  requirement  from  dynamical  considerations 


Prom  the  figure. 

-Mer -Mr  =0  (20) 

where  the  individual  terms  are  moments  about  the  flapping  hinge 
caused  by  the  lift,  blade  weight,  centrifugal,  and  inertia 
forces,  respectively.  Some  of  these  can  be  evaluated  by  in- 
spection: ft 

- / Mgrdr  - mg  ^ <2» 

*r(r/5)dr  < 22 , 

= fo  n(-&0  ’-/4 ) A Vr 

= m(£+ 20.(1*) £% 

+tn££{2*  (23) 
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Mj--  £ (fnr/6)tc/r  = M/$  (24) 

The  lift  moment  is.  derived  from 

f*r  # dr  (25) 

di-^{d-J^)Urdr  (26) 

Substituting  as  before  and  integrating, 

Ml  m -f  [ % 4 (& +2  Q.  Q)  +4 

-GoM  +8,/e<ur-Ct  (27 ) 

Consequently,  equation  (20)  can  be  expanded  to  give,  for  the 
variation  in  flapping  angle  ft, 

* eA6,Ql 

-tUA  v 6,fi+e,»  -c,  yc)  -4,  qj 

- mg. 4 -/n(fc+zQ.  QJfe 

-m/3  e/3  - o (28) 

where  it  might  be  noted  that 

ft  =/i  , and  ft  . 
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IV  VERTICAL  MOTION  OF  THE  HELICOPTER 


"The  equation  for  vertical  motion  of  the  helicopter  is 
simply: 

-jjf  7‘  <291 


w-r 


for  which  only  the  value  of  rotor  'thrust*  is  required  to  permit 
the  solution  to  be  obtained.  In  this  case,  the  rotor  'thrust' 
is  taken  to  be  the  net  force  transferred  to  the  helicopter  by 
the  blade  system.  Thus,  from  figures  (1)  and  (2),  the  elemental 
•thrust'  is  given  by 

</T~ctl-  <pdD  -dl  l30) 

where  the  blade  weight  has  been  omitted  since  it  is  considered 
as  part  of  the  helicopter  weight  W in  equation  (29) • As  is 
well  known,  can  be  neglected  in  evaluating  equation  (30) . 

Writing  c/1  as 

dl  * mrfio/r  - . od 

and  combining  with  equation  (26),  integrating  and  substituting 
in  equation  (29),  there  results 


\o0 C,  (Q0  ■* zQ0Q &) + QjD.0 C, 

-Q> (D,  +C, /S  +F,ur 
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V SUMMARY  OP  EQUATIONS 


The  torque,  flapping,  and  vertical  motion  equations, 
equation  (13),  (28),  and  (32),  respectively,  constitute  a 
set  of  three  simultaneous  differential  equations  and  can  be 
written  in  the  following  form: 


Q/)*Qo  +QurU>-  +Qfi/i  * QeA  * #3*0*  +#jc  f‘ 

0 = Me  * Merit/"  *Ml  p & / MnL  fti 

ft 

0=To  *7^-ur*T^A  * TSa  9a  + £s,ft 
*rjt  i +-7jch  + T/tfi 


(33) 


J ■91,  \ 


(35)  - 


The  first  term  on  the  right  hand  side  of  these  equations  repre- 
sents the  steady-state  term,  and  for  purposes  of  solving  equa- 
tion (1),  each  of  these  can  be  ignored.  However,  these  may 
be  of  interest  in  establishing  the  values  of  steady-state 
flapping  angle,  collective  pitch  angle,  etc.  The  coefficients 
of  the  variables  in  equation  (33)-(35)  are  obtained  from  the 
basic  equations  (13) , (28),  and  (32),  and  together  with  the 
steady-state  terms  are  listed  below: 

Qc  - *(4>&4t-£,)+  %<£&*$&*&€) 

Qv  - a.(0oQ  C, +2Sz0o)QCi  +2%% 
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Q/i  --  a(A Q* 8,-24)- (Si  +2£t Q)JZ 4 +2 6*  A 

(pQ  * aQo  A/  -t-  8 Z (Si  / 2S2  @0  )SO  - 2SzQa  Az 
Qq^~  A Ga  A,  + 2et(£0  +3, 4 +&z  &}Q 
- (S  + 2 £2  Se)  Az 

Q%  - -a  (AM -2 a)  + (6^  2Sz&)Q,C. 

-ZDzSz 

Mo‘  jpf&SQo  '~-Qo  4)  -rtgtfz  -r/iQ ^ S/3  =3 
/%-=  -aQoC./t 

M/S = -^z^>  4-4 
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- a $L I'A 

=•  Z a.&,&eQa/t>  +a.A/t>  -2/n ^ 

/Mjc  - a.Q0  C,/i> 

4fe-*  -/n^}3 

Te  - m/-£4,c,C£'-CLq 

= Q-oh  ^ 

T/t*  QeC,a 
7%  ’-a02C, 

Tn^  * -aQO/ZOo  +al>, 

T^-.-aCLr, 

^ ’ " ?V 
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VI  SYMBOLS 

a slope  of  lift  coefficient  curve,  1/rad. 
b number  of  blades  in  rotor 

B tip  loss  factor 

c blade  chord,  ft. 

CDq  drag  coefficient 

lift  coefficient 

dD  drag  force  on  the  blade  element,  lbs. 
dL  lift  force  on  the  blade  element,  lbs. 
d( CP) centrifugal  force  on  blade  element,  lbs. 
dWb  weight  of  blade  element,  lbs. 
dl  inertial  force  on  blade  element,  lbs. 
dT  blade  element  thrust  (see  equation  30) 
g acceleration  of  gravity,  ft/sec^ 

Iefj.  moment  of  inertia  of  rotating  system,  ft.  lb.  sec^ 

M(jp  centrifugal  force  moment  about  flapping  hinge,  ft.  lb 
nJj  inertia  moment  about  flapping  hinge,  ft.  lb. 
ml  aerodynamic  moment  about  flapping  hinge,  ft.  lb. 

“Wb  blade  weight  moment  about  flapping  hinge,  ft.  lb. 
m mass  density  of  blade  per  foot  of  span,  slugs/ft. 

change  in  rotor  aerodynamic  torque,  ft.  lb. 

Qe  change  in  available  engine  torque,  ft.  lbs. 
r radius  of  a blade  element,  ft. 

R rotor  radius,  ft. 

T integral  of  dT  for  b blades. 
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